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1.  GENERAL  PRESENTATION 

1.1.  Introduction 

The  test  case  RCM-3  consists  of  modeling  the  MASCOTTE  combustor  at  a pressure  chamber  of 
60  bar,  which  is  higher  than  the  critical  pressure  of  liquid  oxygen  (50.4  bar).  This  range  of  pressure  is 
thought  to  be  representative  of  the  chamber  pressure  encountered  in  real  engine. 

The  CPS  code  (version  1.3)  was  used  to  model  this  test  case.  Various  models  have  been  compared 
between  them  and  to  the  experimental  results  available. 


1.2.  Geometry  of  MASCOTTE 

The  MASCOTTE  test  combustor  has  a square  section  of  50  x 50  mm.  The  chamber  length  is  fixed  at 
400  mm. 

The  injector  head  consists  of  a single  coaxial  injector  element  (Figure  1 ). 
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Figure  1 : injector  geometry 


1.3.  Operatins  condition 

The  operating  point  is  the  A-60  case.  It  is  defined  in  the  following  table: 


Pressure 

O/F 

m (LOx) 

m (H2) 

60  bar 

1.4 

100  g/s 

70  g/s 

1 
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In  this  case  [1],  liquid  oxygen  is  injected  at  85  K,  while  gaseous  hydrogen  is  injected  at 
approximately  287  K.  The  physical  properties  of  the  propellants  are  summarized  in  the  following  table: 


Conditions 

H2 

O2 

Pressure 

6 Mpa 

6 MPa 

Mass  flow 

70  g/s 

100  g/s 

Temperature 

287  K 

85  K 

Density 

5.51  kg/m^ 

1 177.8  kg/nf 

Cp 

15110J/kg/K 

1660.9  J/kg/K 

Velocity 

236  m/s 

4.35  m/s 

Surface  tension 

8.67  10"*  kg/m/s 

2.34  10"*  kg/m/s 

1.4.  Experimental  results 

Available  data  for  this  test  case  is  OH  emission.  Abel  transform  permit  to  have  an  axisymetrical 
view  of  the  flame  (Figure  2) 


Average  OH*  emission  image  for  operating  point  A-60 


Abel  transform  emission  image  for  operating  point  A-60 
Figure  2 : Experimental  result  for  operating  point  A-60 
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1.5.  Physical  vhenomena 

The  phenomena  involve  in  such  case  are  summarize  in  the  Figure  3. 


1.6.  Oxygen  thermodynamics  properties 


Various  thermodynamics  properties  of  oxygen  at  60  bar  are  plot  on  Figure  4.  We  observe  important 
variation  between  injection  temperature  (85  K)  and  boiling  temperature  (1 19  K). 

Critical  temperature  and  pressure  are  respectively  154  K and  5.04  MPa. 
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Figure  4 : Evolution  of  oxygen  thermodynamics  properties  between  70  et  250  K at  60  bars. 
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2.  MODELING 


2.1.  Mesh 

The  mesh  is  composed  of  37  x 58  cells  as  represented  in  Figure  5. 


Figure  5 : Mesh  (representation  dilated  by  5 on  vertical  direction) 


2.2.  Conditions 

Calculation  have  been  made  with  constant  parameter  and  model  as: 

- Lagrangian  approach 
Explicit  scheme 

- Droplet  fragmentation  model  (TAB) 

- Droplet  vaporization  model 

- Without  droplet  collision 

- Initial  droplet  diameter  (50  pm) 

Other  parameter  and  model  have  been  compared: 

2.2.1.  Injection  condition 

Injection  of  liquid  in  critical  condition  is  not  well  known.  It's  clear  that  injection  of  one  size  droplet 
is  not  representative.  So  we  have  used  different  injection  condition  for  oxygen: 

- 

- Liquid 

- Mixture  with  Gas  / Liquid  ratio  variable 


2.2.2.  Turbulence  model 

Two  turbulence  model  have  been  used: 

- q-O)  Coakley 

- k-e  Jones-Launder 
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2.2.3.  Combustion  model 

Two  combustion  model  have  been  used: 
Coherent  flame  model  (CFM) 

- Eddy  Break  Up  model  (EBU) 


2.2.4.  Oxygen  thermodynamics  properties 

The  oxygen  thermodynamics  properties  have  take  constant  in  our  calculus.  To  take  into  account  the 
droplet  temperature  increase  we  have  take  different  thermodynamics  properties  for  the  oxygen. 


Density  (kg/m ) 


Viscosity  (kg/m/s) 


Surface  tension  (N/m 


Iniection  (85K) 

Boiling  (119K) 

1177.8 

1004.4 

2.34  10'^ 

1.05  10'^ 

1.32  10'^ 

5.82  10'^ 

987.5 

880.6 

Test  case  RCM-3  usiiu'  CPS 


2.3.  Different  case 

The  different  case  test  are  summarize  in  the  following  table  (the  variable  parameter  are  in  italic); 


Injection 

Combustion  model 

Turbulence 

model 

Oxygen  thermodynamics 
properties 

Gas 

EBU 

q - CO 

Injection  (85  K) 

IW/KSSBSMKM' 

CFM 

q - CO 

Injection 

Mixture 

Mixture  ratio 

CFM 

k-e 

Boiling  (1 19  K) 

Liquid 

CFM 

q - CD 

Injection 

EBU 

EBU/Arrhenius 

Mixture 

EBU 

mm 

Injection 

CFM 

Mixture 

CFM 

q - (0 

Injection 

k-e 

Liquid 

CFM 

q - (0 

Injection 

Boiling 

Cv  / p/  enthalpy  / fi 

Mixture 

CFM 

B 

Injection 

Boiling 

Surface  tension 

3.  RESULTS 

3. 1.  Injection  influence 

We  observe  an  important  variation  of  the  flame  length  with  the  ratio  gas  / liquid  injected  (Figure  6 
and  Figure  7).  The  flame  length  increases  with  the  liquid  fraction  (Figure  6).  This  phenomenon  is  due 
to  the  time  increased  necessary  to  vaporize  all  the  droplets.  We  also  observe  a modification  of  the 
reactive  zone  on  the  near  field  of  the  injector  due  to  the  presence  of  gaseous  oxygen. 

I» 

D.3874E+D5 
D.?6e4EtD5 
D.1B4DZ+D6 
D.3147E+D6 
0.5l08Zi-D6 
^ 0.6D50Zt06 
D.U48ztD7 
* D.lJ08Zf07 
D.2?DlZtD7 
D.«l39DEt07 

I0.6S24Z+07 
0.  W7IZt07 
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Figure  6 : Reactive  rate  (mol/kg/s)  : CFM,  q - CO;  oxygen  (injection) 


If  we  decrease  drastically  the  liquid  fraction,  the  flame  length  increases  (Figure  7).  This  is  due  to  the 
increase  of  the  oxygen  injection  speed  to  obtain  the  correct  amount  of  oxygen  injected. 
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Injection  Gas  - Liquid  (1:1) 

Figure  7 : Reactive  rate  (mol/kg/s)  : CFM,  k - £;  oxygen  (injection) 


We  observe  on  Figure  8 that  liquid  oxygen  volume  fraction  can  be  find  greater  that  one  (red  part  on 
the  figure).  This  is  physically  incorrect  and  due  to  the  fact  that  CPS  neglect  the  volume  of  the  ^oplet. 
To  minimize  this  effect,  in  most  case,  we  have  take  a gas-liquid  mixture  for  oxygen  injection. 
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Figure  8 : Liquid  oxygen  volume  fraction 


3.2.  Turbulence  model  influence 
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Figure  9 : Reactive  rate  (mol/kg/s)  : CFM.  Gas  - Liquid  (1:1),  oxygen  (injection) 
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3,3.  Combustion  model  influence 


We  observe  important  variation  of  the  flame  structure  between  the  two  models  figure  10  and 
Figure  11).  With  Gas+Liquid  injection:  EBU  model  gives  a better  flame  structure  in  the  injector  near 
field.  But  the  CFM  gives  better  results  for: 

- Temperature  value  (Figure  12) 

- Flame  structure 
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Figure  10  : Reactive  rate  (mol/kg/s)  : Gas  - Liquid  (^1:1),  k-  £,  oxygen  (injection) 
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Figure  11 : Reactive  rate  (mol/kg/s)  : Liquid,  q - (O,  oxygen  (injection) 
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Figure  12  : Temperature  radius  profile  (x  = 50  mm)  : Liquid,  q - (O,  oxygen  (injection) 


3,4.  Oxygen  thermodynamics  properties  influence 

When  we  change  oxygen  thermodynamics  properties,  we  don’t  observe  difference  in  the  results 
except  in  the  near  field  of  the  injector.  This  is  due  to  the  modification  of  the  surface  tension  and  the 
increase  of  brealcup  with  decreasing  a (a  = 1.32  10*^  N/m  at  85K  et  5.82  10"^  N/m  at  1 19K). 
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Thermodynamics  properties  at  85  K (Injection) 
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Thermodynamics  properties  at  1 19  K (Boiling) 

Figure  13  : Reactive  rate  (mol/kg/s)  : CFM,  k-  e,  Gas -Liquid  ('^1:1) 


In  supercritical  condition  the  surface  tension  drastically  decrease.  To  observe  the  influence  of  o, 
with  have  made  calculation  with  an  arbitrary  value  for  o witch  is  10"^  N / m.  The  flame  length  decrease 
with  the  surface  tension  (Figure  14). 
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Figure  14  : Reactive  rate  (mol/kg/s)  : Liquid,  q - (O,  oxygen  (injection  except  a) 
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3.5.  Global  case  result 


This  case  correspond  to  the  following  condition; 

- Liquid  injection  for  oxygen 

- Turbulence  model:  q - co 

- Combustion  model:  CFM 
Surface  tension  <j  = 1 10'^  N/m 

- Oxygen  thermodynamics  properties  take  at  85  K (injection) 


Figure  17 : O2  mass  fraction 
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Figure  18  : H2O  mass  fraction 


3. 6.  Comparison  with  experimental  results 

If  we  compare  the  axisymetrical  OH*  emission  image  obtained  experimentally  with  the  most 
comparable  calculated  variable  (the  reactive  rate),  we  obtain  Figure  19.  We  observe  a good  agreement 
in  the  near  field  of  the  injector,  but  the  second  part  of  expansion  of  the  flame  is  not  obtained  with 
calculation. 


0 0.05  0.1  m 


Axisymetrical  average  OH*  emission  image 


D D.5DI1E-D1  D.l 

Calculated  reactive  rate 

Figure  19 : Comparison  experiment  — calculation 
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Experimentally  [2]  observe  an  expansion  angle  of  3.4°  for  OH*  emission  and  5.1°  for  HO  emission. 
If  we  measure  the  reactive  rate  angle  obtained  by  calculation,  we  have  4.0°  figure  20).  This  value  is 
coherent  with  experimental  result. 


Figure  20  : Calculated  reactive  rate  expansion  angle 


4.  CONCLUSION 


The  calculations  performed  for  this  test  case  RCM-3  show  important  difference  on  obtained  results. 
Nevertheless,  they  are  in  good  agreement  with  experimental  results  even  if  they  are  few  for  this 
operating  point  (A-60). 

The  lagrangian  approach  for  oxygen  generates  some  interrogation  when  it's  used  in  supercritical 
condition.  It  appears  that  a real  gas  law  for  oxygen  should  be  interesting  for  eliminates  all  the  problems 
relative  to  droplet  injection. 
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